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Abstract 

Proton data from the GOES 6 and 7 satellites and heavy ion data from the IMP-8 satellite have been 
compared to the expected results of Nymmik’s new model for solar particle event fluences. This 
model calculates the energy spectra of ions for protons through nickel for solar particle events, based 
upon the observed proton integral fluence above 30 MeV. Based upon 27 observed proton events of 
solar cycle 22, and three large historical events, with integral fluences above 30 MeV of greater than 
10 6 particles/cm 2 , a reasonable agreement with model predictions is seen for more than half of the 
events. However, several events show a marked departure from the model predictions, leading to the 
conclusion that there may exist more than a single class of event, or that it may be necessary to 
include additional parameters within the model, such as solar disk position of the source flare, or 
height of disturbance in the solar corona. Data for heavy ions, (oxygen and iron), were limited to a 
total of six solar particle events, of which only two occurred in solar cycle 22. The agreement 
between data and the model predictions appeared to be quite good, however this agreement was 
sensitively dependent upon the value taken for the proton fluence above 30 MeV. 


Introduction 


The Model 

During the past few years, a model to describe solar particle events has been developed by Nymmik at 
Moscow State University, (Nymmik, 1997). This model includes descriptions of the energy and/or rigidity 
spectra of both protons and heavy ions up to nickel. In contrast with earlier work, which suggested an 
exponential distribution, this model proposes a power law in rigidity above 30 MeV/nucleon, and a 
flattening of the spectra below that energy. Apart from the fact that a number of natural phenomena can be 
described by power laws as a result of chaotic or non-linear dynamic behavior, there is no physics per se 
upon which this model is based. The values for coefficients and exponents are based upon observations of 
solar particle events, and are therefore statistical in nature. 

Solar particle events are generally described and classified by the size, or event integrated fluence above 
some energy, depending upon the characteristics of the detector(s), and/or the objectives of the researchers 
human rated spacecraft currently provide shielding against particles with energies below about 30 MeV, so 
it is events with significant numbers of particles above this energy which are of primary importance for this 
study. For solar particle events, this event integrated fluence F is given by: 

oo 

F(>30MeV) = J $ (E) dE (1) 

30 

where E is the kinetic energy, and G> (E) is the differential fluence in protons/cm 2 -MeV. In the proposed 
model , (E) is given by: 

-Y 

$(E)dE = C*(p/p 30 ) * dE/P 


( 2 ) 



where p is the momentum, p, 0 = 239.15 MV is the momentum, (rigidity, p/Z), corresponding to a proton 
energy of 30 MeV. For kinetic energies above 30 MeV, the spectral index Y is calculated in the model as 
a function of event size: 


f 6.1 F(>30) < 10 6 
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y = \ 5.3 - 0.8 * sin[ n/2 * {( log(F(>30)) - 7)}] ltf < F(>30) < 10 8 (3) 

1 4.5 F(>30) > 10 8 


For energies below 30 MeV, the spectra flatten, again as a function of the event size. The spectral “droop” 
index a for the lower energies is given by: 


0.75 

a = 0.2 * Yp * ( Y / 4.5 ) 


(4) 


where Yp = 1.0 for F(>30) >3.0* 10 7 

and Yp = F(>30) / ( 3.0 * 10 ) 7 for F(>30) < 3.0 * 10 7 

The coefficient C is given by: 


F(>30) 

C (5) 
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Similarly, the heavy ion spectra are described by power laws in rigidity. The spectral index Y (z) is 
given by: 

y(z> = ’k*Y (P > (6) 

where Y <p> is the index for the protons, and K is a random value of a log normal distribution with mean: 

log )0 K = 0.1 orlc = 1.26 +/- 0.07 


The coefficients C5 z) = C p w * C (p) are calculated from tabulated values of C p (z) . As with the proton 
spectra, there exists a “droop” index which describes the shape of the heavy ion spectra below 30 
MeV/nucleon. This index is given by: 
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where (A IQ) is the mass-to-charge ratio, and is also tabulated in the literature, (Nymmik, 1997). The above 
expressions are for values of Z between 2 and 28, although as will be seen below, the solar event data for 
most heavy ion species are extremely limited. 


The Data 

The last solar cycle (solar cycle 22) has provided researchers with almost complete satellite coverage of 
solar particle events with fluences above 10 6 particles/cm 2 and particle energies above -1 MeV/nucleon. 
Approximately 27 solar proton events in this solar cycle have been recognized, although there is some 
doubt about the duration of some of these, and whether they represent single or multiple events. The 
integral energy spectra above 1 MeV for 21 of these events has been provided to us kindly by Dr. Herbert 
Sauer of NOAA, (Sauer, 1997). Six addition events are taken from earlier published literature, (Sauer, 
1993; Tylka, et al, 1996). The 21 events provided by Dr. Sauer have integral fluences for energies above: 1, 
5, 10, 30, 50, 60, and 100 MeV as does one event from Tylka, et al. In addition, Sauer, 1993 provides 
integral fluences above 355 and 685 MeV for eight of the 21, as well as for five other events. It is noted 
that the refinement of data between the 1993 Sauer publication and the recent set of 21 events does not lead 
to significant differences, and hence the more recent analyses are used for energies < 100 MeV. The solar 
particle events used in the analyses are listed in Table 1, which also includes the positions of the 
corresponding solar flare if available. 

Data for heavy ions has been taken from Tylka, et al, 1997. This includes spectra for two events in solar 
cycle 22. These data were obtained using the IMP-8, VLET, and Galileo spacecraft. The energies range 
from about 5 MeV/nucleon to about 1-2 GeV/nucleon. Four earlier events from solar cycle 21, are taken 
from Mazur, et al, 1992, and from the Chicago/CRT, and have event integrated fluences from below 1 
MeV/nucleon to -300 MeV/nucleon. In two cases, this energy range is also extended to beyond 1 
GeV/nucleon. These events are listed in Table 2, 


The Analyses 


Results 

Using the values given in Dr. Sauer’s data for F(>30), the event integrated proton fluences F(>E) were 
calculated using Nymmik’s model for kinetic energy in the range 5 MeV to 1000 MeV for each of the 27 
events. In referring to specific events, the notation yymmdd is used throughout to identify the event in 
question. The model curves were then superimposed on the corresponding original data. Because the model 
contains a size dependent index, it was instructive to examine the aggregate set of proton events in three 
fluence bins: small: < 10 7 , medium: 10 7 to 10 8 , and large: > 10 8 . With some notable exceptions, it would 
appear that a flattening of the energy spectra does occur with increasing size, as indicated by the model. In 
fact, for energies of 30 MeV and greater, reasonable to excellent agreement is attained between data and 
model for more than half of the events in the set of 27. Figures 1, 2, 3a, and 3b show the kinetic energy 
spectra of the set of 27 GOES proton events. 


Considering first the set of large events, ( F(>30) > 10 8 protons / cm 2 ), without even examining the 
model results, two events stand out as different. Events 891127 and 910323 would appear to have larger 
spectral indices, i.e. to fall off much more rapidly below 100 MeV than the other seven events in that class. 
By including the two high energy values, 910323 appears to maintain the steep spectrum, whereas 891127 
hardens considerably, thus exhibiting a high energy tail which is unlike any of the other solar proton events 
examined. It is suggested that this tail is due to galactic cosmic ray protons which were not removed from 
the data. Comparing the model results with the data for these two events, it clearly does not provide an 
adequate description of the observed spectra. However, the model predictions for the remaining seven 


events in the class are considerably better, particularly for the huge events 890929 and 891019. Figures 4 
through 7 show the model fits to the data, with the proton spectral index Y +/- la, for events 890929, 
891019, 891 127, and 910323, respectively. 


The medium class of events, with F(>30) integral fluences between 10 7 and 10 8 particles/cm 2 exhibits 
a slightly different grouping of spectral indices. The class contains a total of just six events, two of which, 
910513 and 920625 would appear to have flatter spectra throughout their energy ranges than the other four 
events. In fact, it would be very difficult to distinguish the values of the spectral indices for these two 
events from those applying to the large class events. It is interesting to note however, that these two events 
are the only two in the medium class for which the high energy fluences F(>355) and F(>685) are 
available. At least partially due to these high energy data, the fits of Nymmik’s model appears to be the 
best for these two events. The remaining four events form a rather tight group, for which the spectra appear 
to be steeper than the corresponding group of large events, and to be steeper than the model predicts for 
events of this size. As was seen within the large class events, there appears to be two distinct spectral 
shapes; one which Nymmik’s model describes quite well, and one for which the data fall off much too 
rapidly at higher energies. Figures 8 and 9, showing events 910513 and 900319 illustrate these two types of 
event. 


As would be expected, the “small” class, for which F(>30) is in the 10 6 range, has the most number of 
events, twelve in all. Here too, there is an interesting separation into two subclasses. Four of the events, 
881213, 890520, 890725, and 891115 appear to have much flatter spectra above 10 MeV than even the 
main group of large events. The model therefore fails to fit any of these four events. Figure 10, event 
881213 is representative of this group. Again, it is noted that it is only for these four events that the high 
energy data is available, and also, that they are not included in the Sauer, 1997 compilation. Two other 
events in the small class appear to be fit only marginally by the model: 871 108 and 890501. In both cases, 
the spectrum appears to be too flat at higher energies, but without the data at 355 and 685 MeV it is 
uncertain whether this is in fact the case. The remaining six events in the class are well described by the 
model, as is illustrated in figure 1 1, event 910825. 


In addition to these 27 events, three of the largest historical events from previous solar cycles were 
examined. These are: February 23, 1956, November 12, 1960, and August 7, 1972. Details of these events 
are given in Shea and Smart, 1990, and Wilson, et al, 1997. The spectra are modeled as exponentials in 
Wilson, et al, 1993. It is seen that the event of August, 1972 falls ofT very rapidly with increasing energy, 
and clearly cannot be described by the model, (Figure 12). The February,1956 event was observed to be 
very hard at high energies. It is conceivable that this high energy tail is similar to some of the smaller later 
events, however it still appears to be harder than the model spectral index Y = 4.5, as seen in figure 13. 
The event of November, 1960 appears to be the closest of the three to being described by the Nymmik 
model, (Figure 14). 


To summarize the solar proton data analysis, the spectra of 16 of the 30 events analyzed were fit well by 
the Nymmik model. The remaining 12 events were less well fit, in most cases showing characteristics which 
could not be explained within the context of the model, and which indicate that two or more classes of 
spectral shape may be present. It would seem that additional physical parameters are necessary in order to 
adequately describe the full range of events. It was thought originally that the observed variation in spectral 
shapes might be due to the position on the solar disk of the causal flare. An examination of available data 
showed little or no correlation between flare position and spectral characteristics. Several authors, 
(Reames, et al, 1994; Kahler, et al, 1997), have suggested that temporal variations play a role in the 
composition and the spectral shape, and have noted that source height variation in the corona may also be of 
significance. 


Heavy Ions 


It has been assumed in the past that the heavy ion component of solar particle events was insignificant so 
far as radiation dose calculations were concerned. Additionally, very few measurements have been made at 
satellite altitudes. Early estimates of the composition and spectral characteristics have been done by Mazur, 
et al, 1992, and Reames, et al, 1994. More recently, Tylka, et al, 1996, 1997; have published some spectral 
data for oxygen and iron ions for six solar particle events using the IMP-8 spacecraft, as well as several 
others mentioned above. Following the previous notation, these events are designated: 771122, 780428, 
790819, 810424, 890929, and 891024. It will be noted that two of these events correspond to proton events 
of solar cycle 22. The event 891024 represents a portion of 891019. 

In order to obtain event integrated fluence spectra for the ion spectra, it is necessary to determine the 
corresponding proton fluence F(>30). For the four events in solar cycle 21, this was done using the data in 
Shea and Smart, 1990. Care had to be taken to align the temporal periods of data collection, as given in the 
two sources, in order to get the best estimate for the proton fluence corresponding to the heavy ion event. 
The best illustration of this is found in event 891024. Examining the original data given in the weekly 
NOAA Solar Geophysical Data publication, it was seen that event 891019 actually consisted of three pulses 
of particles. Only the last, (and smallest) pulse coincides with 891024. Therefore, the value of F(>30) for 
891019 had to be corrected. This was done by hand-integrating the temporal flux curve. This integration 
gave a value of F(>30) = 5. 55x1 0 8 for the period October 24-27, compared to 4.24x1 0 9 for the entire 
event. Figures 15a and 15b show the spectra for oxygen and iron particles respectively, fit with the mean 
calculated from the proton fluence F(>30) for this event. A correspondingly good fit was obtained for the 
September 29 th , 1989 event. 

The model results for the two solar cycle 22 events were very encouraging for both oxygen and iron 
spectra. In the published data for these events there is almost no indication of a flattening below 30 
MeV/nucleon. Consequently, we ran the model without the droop calculation for energies below this 
energy. It is observed that the fits are excellent. One possible improvement in the model might be to 
perform the droop calculation with an expression eliminating the A/Q dependence, which seemed to drive 
the low energy portion of the spectrum far too low. 

The model results for the other four events were not quite as good. For 771122, the oxygen model 
paralleled the data about an order of magnitude too high, until intersecting what are possibly galactic 
cosmic ray (GCR) values at approximately 200 MeV/nucleon. The iron spectrum appeared to match the 
data quite reasonably above -20 MeV/nucleon. This would indicate that the ratio of O/Fe is not correct, or 
varies from event to event as indicated in the references. The oxygen data for event 780428 fell off more 
rapidly with increasing energy than the model predicts, even appearing to be somewhat like the exponential- 
type proton spectra. The corresponding iron spectrum has a similar form, although not so extreme. The 
model calculation presents a barely acceptable fit to the data. Both the August 1979 (790819), and the 
April, 1981 (810424) events presented good model fits for both oxygen and iron data. 

To summarize the heavy ion model results, very good fits were obtained in four out of the six cases. The 
model iron results appeared more robust in the remaining cases. It should be noted again that no droop 
calculations were used in the heavy ion spectra model; and that the results are quite sensitive to the values 
taken for the event integral proton fluences. In some cases there also appears to be uncertainty concerning 
whether the data represent GCR or solar particles, particularly at energies above -200 MeV/nucleon. This 
could indicate that the solar heavy ions do indeed represent a more serious radiation dose hazard to 
astronauts than has been assumed previously, particularly in cases such as the October 1989 event, in which 
possibly only part of the dose due to heavy ions was observed. This needs to be examined carefully in the 
near future. 
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Table 1 Proton Events 


Event ID 
(yymmdd) 

Onset Date 

Integral Fluence, F(>30) 
(protons/cm 2 ) 

Solar Disk Position of! 

560223 

February 23, 1956 

1.0 x 10 9 


601112 

November 12, 1960 

9.72 x 10’ 


720804 

August 4, 1972 

5.0 x 10 9 


871108 

November 8 , 1987 

1.20 x 10 6 

N31 W >90 

881213 

December 13, 1988 

8.52 x 10 4 

N27 E37 

890308 

March 8 , 1989 

4.91 x 10 7 

N35 E46 

890411 

April 11, 1989 

4.69 x 10 6 

N37 W12 

890501 

May 1, 1989 

1.28 x 10 6 

S18 W0 

890520 

May 20, 1989 

5.96 x 10 6 

S20 E5 

890725 

July 25, 1989 

8.11 x 10 6 

N28 W84 

890812 

August 12, 1989 

1.52 x 10 9 

S18 W36 

890929 

September 29, 1989 

1.39 x 10 9 

S25 W >90 

891019 

October 19, 1989 

4.24 x 10 9 

S27 E10 

891115 

November 15, 1989 

5.70 x 10 4 

Nil W33 


891127 

November 27, 1989 

1.31 x 10* 

N25 W3 

900319 

March 19, 1990 

2.17 x 10 7 

N33 W43 

900428 

April 28, 1990 

4.80 x 10 6 

unknown 

900521 

May 21, 1990 

1.38x10* 

N35 W36 

900801 

August 1, 1990 

7.41 x 10 6 

N20 E45 

910131 

January 31, 1991 

1.28 xlO 6 

S17 W35 

910323 

March 23, 1991 

1.78 x 10 9 

S26 E28 

910513 

May 13, 1991 

1.82 x 10 7 

S10 W >90 

910604 

June 4, 1991 

7.85 x 10 8 

N30 E70 

910611 

June 11, 1991 

6.27 x 10* 

N31 W20 

910629 

June 29, 1991 

3.00 x 10 7 

N30 E >90 

910825 

August 25, 1991 

3.41 x 10 6 

N24 E77 

920509 

May 9, 1992 

1.32 x 10 7 

S26 E8 

920625 

June 25, 1992 

4.77 x 10 7 

N9 W67 

921102 

November 2, 1992 

2.64 x 10* 

S25 W >90 

940220 

February 20, 1994 

4.13 x 10 6 

Nil W14 


Table 2 Heavy Ion Events 


Event ID Onset Date Proton Integral Fluence F(>30) 

(yymmdd) (protons/cm 2 ) 

771122 November 22, 1977 6.3 x 10 7 

780428 April 28, 1978 2.47 x 10* 

790819 August 19, 1979 9.5 x 10 7 

810424 April 24, 1981 7.0xl0 7 

890929 September 29. 1989 1.39 x 10 9 

891024 October 24, 1989 5.55x 10* 


Figure Captions 


Figure 1: “Large” solar proton events, with integrated event fluences F(>30MeV) > 10* particles/cm 2 . 

Figure 2: “Medium” events, with integrated event fluences 10 7 < F(>30) < 10 s particles/cm 2 . 

Figures 3a, 3b: “Small” events, with 10 6 < F(>30) < 10 7 particles/cm 2 . The class is divided into two in 

order to show the individual events more clearly. 

Figures 4 through 1 1 : Representative individual events (GOES data), from solar cycle 22, together with 

Nymmik’s model spectra. Also shown are the +/- la r values for the calculated y. 

Figure 12: The August 4, 1972 event. Wilson, et al, (1997), describe the available data for the 

differential fluence, (protons/cm 2 -MeV), by the expression: <I>(E) = 3.0 x 10* exp[-(E-30) / 26.5]. This 
expression is integrated, and plotted together with Feynman, et al’s (1993) data, as is the Nymmik model 
spectrum, with +/- la. 

Figure 13: The February 1956 event. The event data are summarized by Wilson, et al (1997) with the 
expression: O(E) = 6.0 x 10 7 exp[-(E - 10) / 25] + 9.4 x 10 5 exp[-(E - 100) / 320] protons/cm 2 -MeV. This 
expression is integrated, and compared to the Nymmik model spectrum for the fluence F(>30MeV), which 
is also shown. 

Figure 14: The November 1960 event. Wilson, et al, (1997) summarize these data with the expression: 
0(E) = 6.3 x 10 s exp[-(E-10) / 12] + 4.9 x 10 6 exp[-(E-100) / 80]. This expression is integrated and 
plotted, together with the corresponding spectrum from Nymmik’s model. 

Figures 15a and 15b: These plots show the Tylka, et al data for the oxygen and iron particles for the 
October 24, 1989 event. Nymmik’s model spectra for these species are calculated from the corresponding 
proton fluence F(>30MeV), for the period 24-27 October. 
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